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Abstract. 

The formulation of statistical physics using light-front quantization, in- 
stead of conventional equal-time boundary conditions, has important advan- 
tages for describing relativistic statistical systems, such as heavy ion colli- 
sions. We develop light-front field theory at finite temperature and density 
with special attention to quantum chromodynamics. First, we construct the 
most general form of the statistical operator allowed by the Poincare algebra. 
In light-front quantization, the Green's functions of a quark in a medium can 
be defined in terms of just 2-component spinors and does not lead to doublers 
in the transverse directions. Since the theory is non-local along the light cone, 
we use causality arguments to construct a solution to the related zero-mode 
problem. A seminal property of light-front Green's functions is that they are 
related to parton densities in coordinate space. Namely, the diagonal and off- 
diagonal parton distributions measured in hard scattering experiments can be 
interpreted as light-front density matrices. 



Dirac's front form of relativistic dynamics 1 has remarkable advantages in 
high energy and nuclear physics. Most appealing is the simplicity of the vacuum 
(the ground state of the free theory is also the ground state of the full theory) 
and the existence of boost-invariant light-cone wavefunctions (see Ref. for 
a review.) This makes light-front quantization a natural candidate for the de- 
scription of systems for which boost invariance is an issue, such as the fireball 
created in a heavy ion collision or the small-x features of a nuclear wavefunc- 
tion. Until now, however, most applications of Dirac's front form refer to the 
case of zero temperature. It is clearly important to exploit the advantages of 
light-front quantization also for thermal field theory. In this talk, we report our 
recent findings, see Ref. |31, where we applied front-form dynamics to statistical 
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physics and investigated the prospects and chahenges of this approach for quan- 
tum chromodynamic systems. Valuable work in this direction has already been 
done by several authors |^ El El • 

The form of the statistical operator w at finite temperature and density can 
be obtained from very general considerations. Our result for w, is compatible 
with the findings of Ej, id is always the exponential of the equal time 
energy in the local rest frame of the system. Our derivation follows Ref. 0. 

The light-front Liouville theorem {id~w = [P~,w]) requires that in equi- 
librium, w IS a function of only those Poincare generators which commute with 
the light-front Hamiltonian P~ . In addition, since systems far apart from each 
other must be uncorrelated, the density operator of the combined system has to 
factorize into the density operators of the subsystems. Consequently, in equilib- 
rium ln{id) must be a linear combination of the additive constants of motion, 
namely the four components of the momentum and the 3-component of the 
angular momentum vector J. Hence, 

ln{id) = a-(3(u,P'' -uJs-J^l^iQi^ ■ (1) 

Here, f3 is the inverse temperature and Ui, is the four velocity of the system, 
cf. Ref. [H]. In addition, to is the angular velocity at which the body rotates. 
Additional conserved charges Qi are included along with their chemical potentials 

In quantum mechanics, of course, one can simultaneously specify only charges 
which commute with each other, e.g. one cannot specify all four P^ for systems 
with nonzero angular momentum. 

We remark that T is the same as the instant-form temperature, but the 
chemical potential has a different meaning. On the light-front, densities are given 
by -|— components of currents, and not by 0-components. This is essential for 
the proper generalization of parton distributions (PDFs) to finite temperature, 
since the latter are also defined as -|— components. Finite temperature PDFs are 
useful for parton recombination models (see e.g. 0), even though they cannot 
be measured in deep inelastic scattering. 

We choose a = as normalization in Eq. so that the partition function 
is given hy Z = Tiw. The grand-canonical ensemble can now be written in terms 
of light-cone wavefunctions (p^/hi^) &s (let P± = 0_l), 

h n,n' X,X' I 

X 0,/,(XX,/,(X')|nX)(n'X'|. (2) 

Since the wavefunctions and the masses Mh of the eigenstates can (in principle) 
be obtained from discretized light-cone quantization (DLCQ) |9j , this expression 
shows that one can also calculate all thermodynamic properties of a field theory 
from DLCQ. 

Furthermore, since Z is a Lorentz scalar, all thermodynamic potentials and 
the entropy transform as scalars, e.g. the Lorentz invariant generalization of 
the grand-canonical potential (or of the free energy in the case of /i = 0) is 



J. Raufeisen and S.J. Brodsky 



3 



1? = —TlnZ{V,T,n), and the entropy is given by 

The role of the total energy of the system is now played by the expectation value 
of UyP^ , U = {uyP^). As usual, Q = U — TS — fiQ. All known relations between 
thermodynamic potentials remain valid. 

The quantities /3, u^, lo and /i/, have the meaning of Lagrange multipliers 
that hold the mean values of the constants of motion fixed, while entropy is 
maximized. In an ideal gas for example, the maximum entropy is attained for 
occupation numbers given by Fermi-Dirac and Bose-Einstein statistics (U 1^1 , 

niu^p") = ^ ^ - g/3(„,p-^+;,) ^ ' (4) 

assuming that particles carry charge +1 and antiparticles charge —1. (Note that 
Uup'^ > 0.) Here, p'^ is the 4- momentum of a single particle and the degeneracy 
factor for different spin states is denoted by g. The Lagrange multipliers define 
the equilibrium conditions for two systems. In complete equilibrium with each 
other, both systems must have the same values of temperature, Uy, uj and fii, i.e. 
no internal motion of macroscopic parts of the system is possible in equilibrium 
(at least in the absence of vortex lines (2|-) 

The simplicity of the light-front vacuum, usually considered an advantage, 
seems to bear problems as far as phase transitions are concerned. However, the 
statistical weight of a configuration is maximized for minimal equal-time energy 
rather than for minimal light-front energy. Therefore, the ground state, i.e. the 
state the system is in at T = 0, is in general different from the light-front 
vacuum. For that reason, the authors of Ref. obtain the standard pattern 
of spontaneous symmetry breaking in (/>^-theory with negative mass squared. 
No problem arises from l/A:"^-poles. In addition, in Ref. |4j the chiral phase 
transition in the Nambu-Jona-Lasinio Model on the light-front is reproduced. 
We conclude that this approach is poised for the study of phase transitions in 
more complicated field theories, such as QCD. 

Until now one could get the impression that thermodynamics and statistical 
physics on the light front are identical to the usual instant form approach, except 
for a trivial change of variables. That this is not the case becomes most clear, 
when one studies fermions on the light-front. In light-front field theory, the Dirac 
equations can be written as a set of two coupled equations for 2-component 
spinors, see e.g. appendix of Ref. Only one of these equations contains a time 
derivative, the other one is a constraint. As a consequence, the entire theory can 
be formulated in terms of 2-component spinors, very much like a non-relativistic 
theory. 

The time-ordered Green's functions of a fermion in a medium is defined in 
terms of the Heisenberg operators of the dynamical field components |K1 llOj. 

iG«,/3(ri,r2) = {Mri)i'l{r2))0{rt - 4) - {i'l{r2)Mri))0{4 - rt), (5) 

where the average (. . .) is to be taken with the appropriate ensemble, and 
a,/3 E {1,2}. This definition of the Green's function includes the case of zero 
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temperature. Therefore, the conventional hght-front quantization at temperature 
T = can be formulated in terms of Ga^p as well. The Green's function is the 
fundamental object of this approach. In addition, the retarded {R) and advanced 
{A) Green's functions are defined by the anticommutators 

i<'^(^i,^2) = ±({^„(ri),V^J(r2)})0(±(r+-r2+)), (6) 

where the upper sign refers to G^ ^ and the lower sign to G^^. We remark, that 
in a gauge theory, it is also necessary to include a (path ordered) gauge link 
along the light-cone by redefining the fermion fields, see Ref. 0. 

The Green's functions of a fermion in an ideal gas of temperature T = 1/(3 
were presented first in Ref. [3]. In momentum space, adjusted to our notation, 
they read 




"'^ " k'^ — m? ± iesgn(uA;) ' 




^SJjW = 5„,^(^P^2^-isgn(^^^Ktanh(^^)A;+5(A;2- 

where P refers to principle value prescription. The pole prescriptions ibiesgn(n/c) 
for the retarded and advanced Green's functions are another manifestation of 
the special meaning of the equal-time energy. These prescriptions ensure that 
G^a^i'^) vanishes outside the forward lightcone, while G^^'^^{r) is non- vanishing 
only inside the backward lightcone. Most importantly, knowledge of the correct 
pole prescription eliminates ambiguities in the definition of the non-local operator 
l/Zc^, which appears in the free light-cone Hamiltonian. However, the correct 
prescription for the l/A;"'"-pole depends on the type of Green's function and on 
the value of the other momentum components. 

Another remarkable property of the light-front Green's functions is, that if 
the theory is discretized on a lattice in coordinate space, the factor k'^ in the 
numerator leads to only one pair of fermion doublers. For a transverse lattice 
lattice approach to finite temperature SU (oo), see Ref. 11 . Moreover, it is known 
that no fermion doubling problem occurs in DLCQ and one can perform DLCQ 
calculations for a fixed value of the total charge without any sign problem. 

In the limit r"*" — > 0^, the time-ordered Green's function Gq^/j is closely 
related to the one-particle density matrices for fermions and antifermions, 
Qa,l3{Li-,L2) ^"^^ Qa,i3iLiiL2)- the 2-component theory, the separation of fermion 
and antifermion distributions requires the evaluation of a Fourier-integral of the 
Green's function. For quarks, one has 

QaAk^^&r^) = -^J dr-e+^^^^'/^G^^fsirt ^ 0-,r„rt,r^). (9) 

For antiquarks, the limit r+ ^ is taken from the other side to obtain the correct 
order of creation and annihilation operators jnj- Since Ga,f3{ri,r2) often depends 
only on the difference r = ri — r2, we introduce the variables R= {ri + r2)/2 and 
r = ri — r2. The density matrix Qq,^^ is related to the so-called Wigner function 
by a Fourier transform over r±. We remark that all properties of the quantum 
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mechanical density matrix, such as hermiticity and positivity of the diagonal 
matrix elements, also apply to the light-front density matrix in = gauge, 
since this gauge has only states with positive norm and no unphysical degrees of 
freedom. However, the light-front density matrix has matrix elements that are 
off-diagonal in Fock-space. The object defined in Eq. Q is similar to the Wigner 
function introduced in Ref. |12j . 

The fermion density matrix contains all information about single-quark prop- 
erties. It depends on 6 variables and is a 2 x 2 matrix in spinor space, which can 
be written as a linear combination of Pauli spin matrices. The coefficients are 
the density matrices for unpolarized, longitudinal, and transverse spin distribu- 
tions. The diagonal matrix elements in coordinate space of (?a,/3(A:''',i?, J'x), i-e. 
the ones with r± = Oj_, are closely related to the usual PDFs. For instance, the 
unpolarized collinear quark density is given by 

q{k+) = ^J cfiR6a,f3qisAk^^&r± = ^±)- (10) 

This parton density is normalized such that /q°° dk^q{k^) = q, the total number 
of quarks in the system. 

The off-diagonal matrix elements of qa^/^ are related to generalized parton 
distributions (GPDs) by Fourier transform The precise relation depends 
on the kinematics, and one has to distinguish four different domains in deeply 
virtual Compton scattering (DVCS). In particular, for skewedness C = 0) GPDs 
can be identified as impact parameter dependent parton densities ^3]. In the 
case of no helicity flip we find, 

q{k+,b) = I dR-hfi,a.qc,,p{k+,R-,h,r^=Ox) (H) 
Y}\k+,h,X)h{k+,h,\)Y (12) 

q{k+,b)dk+ = / r_^e'''-^^//(X,C = 0,t), (13) 

We use the notation of Radyushkin here [T3j . The creation operator of a quark 
at impact parameter h is given by 

b{k+,hA) = J ^e-"'^-^b{k+,k^,\). (14) 

The destruction operator 6^(A;+,6, A) is defined analogously. We stress that for 
C 7^ 0, GPDs are in general not probability distributions but density matrices, 
which do not need to be positive. 

The light-front density matrix is a natural extension of the parton model to 
quantum mechanics: classical parton densities are replaced by a density matrix. 
It would be interesting to identify other hard processes besides DVCS which 
are sensitive to the quantum mechanical nature of parton distributions. Most 
important however is the connection between the density matrix and the fermion 
Green's function, because that establishes a common language for high energy 
scattering and statistical QCD. 
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In summary, we have presented a new formalism for analyzing relativistic 
statistical systems based on light-front quantization. The new formalism pro- 
vides a boost-invariant generalization of thermodynamics, and thus it has direct 
applicability to the QCD analysis of heavy ion collisions and other systems of 
relativistic particles. 
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